dat_general=[];str=[];xs=[];Disp=[];EX_general=[];DX_general=[];
dat_ori=[];
n=35;
v=1;
v=3.644551e-04
dt=1;
shift=200;

for i=5:10:n
 
Pic_Num=5000*i;

if (i==0)
    Pic_Num=30000;
end


    dat=load(strcat('Statistical_data_concentration_X_',int2str(Pic_Num),'.sta'));
    [nx,ny]=size(dat);

dat_ori=[dat_ori,dat];
dat1=[];

    for j=shift+1:nx
        dat1=[dat1;dat(j)];
    end
    for j=1:shift
        dat1=[dat1;dat(j)];
    end
    
dat_general=[dat_general,dat1*v*dt*Pic_Num];
x=1:nx;
x=x-(nx-shift);
x=x/(dt*v*Pic_Num);
plot(x,dat1*v*dt*Pic_Num);
hold on 
%dat_general=[dat_general,dat1];
end

% x=1:nx;
% x=x-(nx-shift);
% plot(x,dat_general);

 
 



